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Abstract. With the goal of confirming the metallicity "excess" present in stars with planetary-mass companions, 
we present in this paper a high-precision spectroscopic study of a sample of dwarfs included in the CORALIE 
extrasolar planet survey. The targets were chosen according to the basic criteria that 1) they formed part of a 
limited volume and 2) they did not present the signature of a planetary host companion. A few stars with planets 
were also observed and analysed; namely HD6434, HD 13445 (G186), HD 16141, HD 17051 (t Hor), HD 19994, 
HD 22049 (e Eri), HD 28185, HD 38529, HD 52265, HD 190228, HD 210277 and HD 217107. For some of these 
objects there had been no previous spectroscopic studies. 

The spectroscopic analysis was done using the same technique as in previous work on the metallicity of stars with 
planets, thereby permitting a direct comparison of the results. The work described in this paper thus represents the 
first uniform and unbiased comparison between stars with and without planetary-mass companions in a volume- 
limited sample. The results show that 1) stars with planets are significantly metal-rich, and 2) that the source 
of the metallicity is most probably "primordial". The results presented here may impose serious constraints on 
planetary system formation and evolution models. 
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1. Introduction 

One of the most exciting and promising results that be- 
came evident after the discovery of the first extrasolar 
planets is that their host stars seem to be very metal- 
rich when compared with dwarf stars i n the s olar neigh- 
bourhood (Gonzalez 1998; Santos et al. P000a|- hereafter 



Paper I; Gonzalez et al. 2001 ). This result, representing 
the only link between the presence of planets and a stellar 
photospheric feature, has been given two main explana- 
tions. The first, is based on the classical view that giant 
planets are formed by runaway accretion of gas on to a 
"planetesimal" having up to 10 Earth masses. In such a 
case, we can expect that the higher the proportion is of 
dust to gas in the primordial cloud (i.e. metals), and con- 
sequently in the resulting proto-planetary disc, the more 
rapidly and easily may planetesimals, and subsequently 
the now observed giant planets be built. 



Send offprint requests to: Nuno C. Santos, e-mail; 
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* Based on observations collected at the La Silla Observatory, 
ESQ (Chile), with the CORALIE spectrograph at the 1.2- 
m Euler Swiss telescope, and with the FEROS spectrograph 
at the 1.52-m ESO telescope (Observing run 66.C-0116B), 
and using the UES spectrograph at the 4-m William Hershel 
Telescope (WHT) , at La Palma (Canary Islands) . 



Opposing to this view, it had been proposed that the 
observed metallicity "excess" may be related to the "pol- 
lution" of the convective envelope of the star by the in- 
fall of planets and/or planetesimals (e.g., Gonzalez 1998; 
Laughlin & Adams 1997; Laughlin pOOO ; Gonzalez et al. 
pOOl| ). This pollution can be the result of the "complete" 



inward migration of a planet on to the star, the transfer 
of material from the disc to the star as a res ult of the 
migration process (Goldreich & Tremaine |l980t Lin et al. 
1996 ), or to the break-up and infall of a planet (s) on to 
the surface of the star due to gravitational interactions 
with other companions (Rasio & Ford |l996D . The former 
point might be particularly important for the short-period 
systems (Queloz et al. 



2000D 



The idea that planets and/or planetary material might 
be engulfed by a star is to some extent supported by the 
recent detection of ^Li in the planet host star HD 82943 
(Israelian et al. 2001) - interpreted as a signal of the accre- 
tion of a planet during the history of the star - and prob- 
ably by the detection of a significant difference in iron 
abundances in the very similar pair of dwarfs 16 CygA 
and 16 CygB, which also have very different Li contents 



(Laws & Gonzalez |2001| ; Gonzalez |1998D . 

Besides the [Fe/H] differences, there is currently some 
debate about possible anomalies concerning other ele- 
ments (Paper I; Gonzalez et al. 2001 ; Smith et al. 2001 ). 
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Table 1. Volume-limited sample of stars without detected giant planets. 



Star 


Tcti 
(K) 


log 3 
(cgs) 


6 

(kms-i) 


[Fe/H] 


N(Fei) 


N(Fell) 


(j(Fel) 


CT(Feii) 


Mass 
[Mq] 


HD1581 


5940 


4.44 


1.13 


-0.15 


31 


7 


0.04 


0.05 


0.99±0.02 


HD 4391 


5955 


4.85 


1.22 


0.01 


36 


5 


0.05 


0.09 


1.22±0.04 


HD5133 


5015 


4.82 


0.92 


-0.08 


36 


6 


0.05 


0.07 


0.81±0.04 


HD 7570 


6135 


4.42 


1.46 


0.17 


35 


7 


0.04 


0.06 


1.19±0.02 


HD 10360 


5045 


4.77 


0.89 


-0.19 


36 


5 


0.04 


0.04 


0.70±0.03 


HD 10647 


6130 


4.45 


1.31 


-0.03 


34 


7 


0.03 


0.05 


1.14±0.03 


HD 10700 


5370 


4.70 


1.01 


-0.50 


38 


6 


0.04 


0.04 


0.68±0.02 


HD 14412 


5410 


4.70 


1.01 


-0.44 


35 


6 


0.04 


0.01 


0.78±0.05 


HD 17925 


5220 


4.60 


1.44 


0.08 


35 


6 


0.07 


0.04 


0.92±0.06 


HD 20010 


6240 


4.27 


2.23 


-0.20 


33 


6 


0.05 


0.08 


1.33±0.01 


HD 20766 


5770 


4.68 


1.24 


-0.20 


35 


7 


0.04 


0.04 


0.97±0.04 


HD 20794 


5465 


4.62 


1.04 


-0.36 


39 


7 


0.05 


0.05 


0.74±0.02 


HD 20807 


5865 


4.59 


1.28 


-0.22 


37 


7 


0.04 


0.04 


0.95±0.02 


HD 23249 


5135 


4.00 


1.12 


0.17 


36 


7 


0.06 


0.08 


0.84±0.01 


HD 23356 


5035 


4.73 


0.96 


-0.05 


36 


6 


0.06 


0.07 


0.83±0.03 


HD 23484 


5230 


4.62 


1.13 


0.10 


37 


6 


0.05 


0.07 


0.92±0.06 


HD 26965A 


5185 


4.73 


0.75 


-0.26 


37 


5 


0.05 


0.03 


0.71±0.02 


HD 30495 


5880 


4.67 


1.29 


0.03 


37 


7 


0.04 


0.03 


1.11±0.04 


HD 36435 


5510 


4.78 


1.15 


0.03 


37 


6 


0.06 


0.03 


1.04±0.05 


HD 38858 


5750 


4.56 


1.22 


-0.22 


36 


7 


0.04 


0.02 


0.91±0.02 


HD 39091 


5995 


4.48 


1.30 


0.09 


37 


7 


0.04 


0.04 


1.10±0.02 


HD 40307 


4925 


4.57 


0.79 


-0.25 


37 


4 


0.06 


0.10 


0.76±0.07 


HD 43162 


5630 


4.57 


1.36 


-0.02 


35 


7 


0.05 


0.04 


0.99±0.04 


HD 43834 


5620 


4.56 


1.10 


0.12 


38 


7 


0.05 


0.05 


0.96±0.02 


HD 50281A 


4790 


4.75 


0.85 


0.07 


31 


4 


0.05 


0.12 


0.79±0.01 


HD 53705 


5810 


4.40 


1.18 


-0.19 


36 


7 


0.03 


0.03 


0.92±0.01 


HD 53706 


5315 


4.50 


0.90 


-0.22 


36 


7 


0.05 


0.08 


0.88±0.05 


HD 65907A 


5940 


4.56 


1.19 


-0.29 


39 


7 


0.04 


0.04 


0.94±0.01 


HD 69830 


5455 


4.56 


0.98 


0.00 


38 


7 


0.05 


0.03 


0.89±0.06 


HD 72673 


5290 


4.68 


0.81 


-0.33 


38 


6 


0.04 


0.04 


0.78±0.06 


HD 74576 


5080 


4.86 


1.20 


0.04 


36 


5 


0.06 


0.05 


0.91±0.04 


HD 76151 


5825 


4.62 


1.08 


0.15 


38 


7 


0.03 


0.04 


1.09±0.04 


HD 84117 


6140 


4.35 


1.38 


-0.04 


34 


7 


0.04 


0.06 


1.14±0.01 


HD 189567 


5750 


4.57 


1.21 


-0.23 


37 


7 


0.04 


0.05 


0.89±0.01 


HD 191408A 


5025 


4.62 


0.74 


-0.51 


37 


4 


0.06 


0.09 


0.61±0.04 


HD 192310 


5125 


4.63 


0.88 


0.05 


36 


6 


0.06 


0.08 


0.78±0.04 


HD 196761 


5460 


4.62 


1.00 


-0.27 


38 


7 


0.05 


0.05 


0.81±0.03 


HD 207129 


5910 


4.53 


1.21 


-0.01 


36 


7 


0.04 


0.03 


1.03±0.02 


HD 209100 


4700 


4.68 


0.60 


0.01 


34 


3 


0.07 


0.06 


0.64±0.09 


HD 211415 


5925 


4.65 


1.27 


-0.16 


35 


7 


0.03 


0.04 


0.99±0.02 


HD 216803 


4647 


4.88 


0.90 


0.07 


28 


3 


0.07 


0.08 


n 77+0.01 


HD 222237 


4770 


4.79 


0.35 


-0.22 


37 


3 


0.08 


0.08 




HD 222335 


5310 


4.64 


0.97 


-0.10 


33 


5 


0.05 


0.04 


0.86±0.05 



But the relatively low number of exoplanets known, and 
possible systematics with respect to the samples do not 
permit firm conclusions to be reached on the subject. 

The observation that stars with planets]^ are particu- 
larly metal-rich has so far been overshadowed by the re- 
striction that in order to compare the metallicities of stars 



^ Here we are referring to the known extrasolar planetary 
host stars, having Jupiter-like planets in relatively short period 
orbits when compared to the giant planets in our own Solar 
System. 



hosting planets with those of stars "without" planets, au- 
thors had access only to published metallicity studies of 
volume-limited samples of dwarfs in the solar neighbour- 
hood (mainly that of Favata et al. (1997), the alternative 
being to construct a sample using less precise metallici- 
ties using photometric indices). This is inconvenient for a 
number of reasons. First and most obviously, one cannot 
be sure if the sample used for comparison is really free 
from giant planets. Furthermore, the metallicities for the 
Favata et al. sample were determined using a much shorter 
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Table 2. Determined atmospheric parameters and [Fe/H] for a set of stars with planets and brown dwarf companions 
(bd). 



Star 


Teff 

(k) 




logs 
(cgs) 




(kms 


') 


[Fc/H] 


Instrument 




Mass 
[Mq] 


HD 1237 


5555±50 


4 


65±0 


15 


1 


50±0 


08 


0.11±0.08 


CORALIE 


1 


01±0 


04 


HD 6434 


5790±40 


4 


56±0 


20 


1 


40±0 


10 


-0.55±0.07 


FEROS 





79±0 


01 


HD 13445 


5190±40 


4 


71±0 


10 





78±0 


10 


-0.20±0.06 


CORALIE 








HD 13445 


5220±40 


4 


70±0 


10 





85±0 


07 


-0.19±0.06 


FEROS 








HD 13445 (avg) 


5205 




4.70 






0.82 




-0.20 







75±0 


04 


tlU lDl41 


00U0±4U 


4 


2o±U 


10 


1 


37±0 


U9 


U.lo±O.Uo 


r liKUo 


1 


06±0 


01 


HD 17051 


6225±50 


4 


65±0 


15 


1 


20±0 


08 


0.25±0.06 


FEROS 


1 


26±0 


02 


HD 19994 


6165±40 


4 


13±0 


20 


1 


49±0 


10 


0.23±0.06 


CORALIE 








HD 19994 


6250±40 


4 


27±0 


10 


1 


56±0 


08 


0.30±0.06 


FEROS 








HD 19994 (avg) 


6210 




4.20 






1.52 




0.26 


- 


1 


34±0 


01 


HD 22049 


5135±40 


4 


70±0 


10 


1 


14±0 


07 


-0.07±0.06 


CORALIE 





85±0 


04 


HD 28185 


5705±40 


4 


59±0 


10 


1 


09±0 


06 


0.24±0.05 


CORALIE 





98±0 


07 


HD 38529 


5675±40 


4 


01±0 


15 


1 


39±0 


09 


0.39±0.06 


FEROS 


1 


52±0 


05 


HD 52265 


6075±40 


4 


21±0 


10 


1 


31±0 


07 


0.22±0.07 


CORALIE 








HD 52265 


6120±50 


4 


37±0 


20 


1 


31±0 


06 


0.25±0.06 


FEROS 








HD 52265 (avg) 


6100 




4.29 






1.31 




0.24 


- 


1 


18±0 


01 


HD 75289 


6135±40 


4 


43±0 


20 


1 


50±0 


07 


0.27±0.06 


CORALIE 


1 


20±0 


02 


HD 82943 


6025±40 


4 


54±0 


10 


1 


10±0 


07 


0.33±0.06 


CORALIE 


1 


15±0 


05 


HD 83443 


5500±60 


4 


50±0 


20 


1 


12±0 


09 


0.39±0.09 


CORALIE 





90±0 


05 


HD 108147 


6265±40 


4 


59±0 


15 


1 


40±0 


08 


0.20±0.06 


CORALIE 


1 


27±0 


02 


HD 121504 


6090±40 


4 


73±0 


10 


1 


35±0 


08 


0.17±0.06 


CORALIE 


1 


19±0 


03 


HD 162020'"* 


4830±80 


4 


76±0 


25 





72±0 


12 


O.OliO.ll 


CORALIE 








HD 168746 


5610±30 


4 


50±0 


15 


1 


02±0 


08 


-0.06±0.05 


CORALIE 





88±0 


01 


HD 169830 


6300±30 


4 


04±0 


20 


1 


37±0 


07 


0.22±0.05 


CORALIE 


1 


42±0 


01 


HD 190228 


5360±40 


4 


02±0 


10 


1 


12±0 


08 


-0.24±0.06 


WHT 





84±0 


01 


HD 202206*"* 


5765±40 


4 


75±0 


20 





99±0 


09 


0.37±0.07 


CORALIE 


1 


02±0 


02 


HD 210277 


5575±30 


4 


44±0 


10 


1 


12±0 


08 


0.23±0.05 


FEROS 





94±0 


01 


HD 217107 


5660±50 


4 


43±0 


10 


1 


04±0 


07 


0.40±0.06 


CORALIE 








HD 217107 


5655±40 


4 


42±0 


05 


1 


11±0 


08 


0.38±0.05 


FEROS 








HD 217107(avg) 


5660 




4.42 






1.01 




0.39 







97±0 


05 



Planet 

discovery paper 



Naef et al. ( pOOlaD 
Queloz et al. ( |2001| ) 



Queloz et al. (200C) 
Marcy et al. (|200C|) 
Kiirster et al. ( |2000[ ) 



Queloz et al. (|200l| ) 
Hatzes et al. (2000) 
Mayor et al. ( 2001b[) 
Fisher et al. ( 2000D 



Naef et al. (gOOlaj) 
Butler et al. (|2000| ) 
Udry et al. (|2000|) 
Naef et al. ( ^OOlaD 
Mayor et al. ( |2001aD 
Pepe et al. (in prep . 
Queloz et al. (|200l|) 



Udry et al. (in prep 
Pepe et al. (in pre p. 
Naef et al. ( bo01a|) 
Sivan et al. ( |2001| ) 
Udry et al. (i n pre p 
Marcy et al. (|l999|) 



Fischer et al. (1999) 



hne list for iron and different sources of atmospheric pa- 
rameters (spectroscopic vs. colours) - see Paper I. This 
particular point may introduce systematic errors, and one 
might expect that the difference between the two samples 
was simply due to a bias related to the method used to 
compute the metallicities. 

With the goal of settling the question about the high 
metallicity content of stars with planets, we present here a 
spectroscopic study of a volume-limited sample of 43 stars 



included in the CORALIE (Udry et al. pOOOj ) planet search 
programme, and for which the radial velocities seem to be 
constant over a large time interval. The technique used, 
line lists and atmospheric models were those usually ap- 
plied by most authors working on the metallicities of stars 
with planets (e.g. Paper I; Gonzalez et al. |2001| ). We show 
that the currently known stars with giant planets are on 
average more metal-rich than "field stars" , for which there 
is no radial- velocity signature of planets. Furthermore, the 
results are used to set strong constraints on the cause of 
the observed "anomaly" , significantly excluding the "pol- 



lution" scenario. The impact on the planetary formation 
and evolutionary models is discussed. 

2. The data 

2.1. The samples 

As mentioned above, the choice of a "comparison sam- 
ple" (as we shall call it from now on), free from giant 
planets and with the minimum bias possible, is of crucial 
importance. In order to be sure that we are not includ- 
ing any systematics when selecting the sample, we have 
taken all stars within the CORALIE sample (defined as 
volume-limited - Udry et al. 20001 ) having right ascen- 
sions between 20*^ and 9*^ (observable in mid November 
from La Silla). Since such a sample corresponds to more 
than half of the ~ 1600 stars searched for planets, it is 
virtually impossible to make a precise spectroscopic anal- 
ysis for all those stars in a reasonable amount time. We 
thus decided as a first step to limit our observations to all 
stars within 17pc of the Sun having {B — V) < 1.1 (deter- 
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Table 3. Stars with giant planets or brown dwarf companions (bd) studied by other authors. 



Star 



(k) 



logs 
(cgs) 



(kms-1) 



[Fe/H] 



Mass 

[Mq] 



Reference 

for spectroscopy 



Planet 

discovery paper 



HD 9826 


6140 


4 


12 


1.35 


0.12 


HD 10697 


5605 


3 


96 


0.95 


0.16 


HD 12661 


5714 


4 


45 


0.99 


0.35 


HD 37124 


5532 


4 


56 


0.85 


-0.41 


HD 46375 


5250 


4 


44 


0.80 


0.21 


HD 75732A 


5250 


4 


40 


0.80 


0.45 


HD 80606 


5645 


4 


50 


0.81 


0.43 


HD 89744 


6338 


4 


17 


1.55 


0.30 


HD 92788 


5775 


4 


45 


1.00 


0.31 


HD 95128 


5800 


4 


25 


1.0 


0.01 


HD 114762'"' 


5950 


4 


45 


1.0 


-0.60 


HD 117176 


5500 


3 


90 


1.0 


-0.03 


HD 120136 


6420 


4 


18 


1.25 


0.32 


HD 130322 


5410 


4 


47 


0.95 


0.05 


HD 1 34Q87 


5715 


4 


33 


1.00 


0.32 


HD 143761 


5750 


4 


10 


1.2 


-0.29 


HD 145675 


5300 


4 


27 


0.80 


0.50 


HD 168443 


5555 


4 


10 


0.90 


0.10 


HD 177830 


4818 


3 


32 


0.97 


0.36 


HD 186427 


5685 


4 


26 


0.80 


0.07 


HD 187123 


5830 


4 


40 


1.00 


0.16 


HD 192263 


4964 


4 


49 


0.95 


-0.03 


HD 209458 


6063 


4 


38 


1.02 


0.04 


HD 217014 


5795 


4 


41 


1.05 


0.21 


HD 222582 


5735 


4 


26 


0.95 


0.02 



1.28±0.02 
l.lOiO.Ol 
1.01±0.02 



1.05±0.03 

~1.1 
1.55±0.03 
1.05±0.02 

1.03±0.03 
0.82±0.03 
1.10±0.02 
1.34±0.02 
~0.89 

1.02-0.03 

0.96±0.03 

~1.05 
1.01±0.02 



1.03±0.11 
0.97±0.03 
1.08±0.04 
~0.80 



1.07±0.01 
0.95±0.01 



Gonzalez & Law s (|200 i 
Gonzalez et al. (2001) 
Gonzalez et al. (2001) 
Gonzalez et al. (2001) 

Gonzalez et al. (2001)^ 

Gonzalez & Vantu re ( |l998| ) 
Naef et al. ( |2001h|) 
Gonzalez et al. (2001) 
Gonzalez et al. (2001) 



Gonzalez 
Gonzalez 
Gonzalez 
Gonzalez 
Gonzalez 
Gonzalez 
Gonzalez 
Gonzalez 
Gonzalez 



( 1998 ) 
1998) 



(|1998D 



Law s (good ) 
et al. ([200 1|) 
et al. (2001) 



( 11998D 

et al. (| l999|) 
et al. (2001) 



Gonzalez et al. (2001) 
Laws & Gonzaler^ra) 
Gonzalez et al. (|1999|) 
Gonzalez et al. (2001) 



1.12±0.02 Gonzalez et al. (2001) 



Gonzalez et al. (2001) 
Gonzalez et al. (2001) 



Butler et al 
Vogt et al. I 
Fisher et al, 
Vogt et al. (I 
Marcy et al, 
Butler et al, 
Naef et al 




1997 



( |2001b|l 



Korzennik et al. (\200(\ ) 
Queloz et al. fcOOlT 
Fisher et al. (GOO 



Butler & Marc y ([1996D 
Latham et al. (|l 98E|) 



Marcy & But ler]l99e| ) 
Butler et al. (|l99 
Udry et al. feOOO l 
Vogt et al. ( pOOl 
Noyes et al. "^99 
Udry et al ( pOOl| 
Marcy et al. (|199 9| 
Udry et al. (|200 l| 
Vogt et al. (pool 



Cochran et a l. ([19971) 
Butler et al (F^ 
Santos et al (jgO OObl 
Vogt et al. ( pOoi 
Mazeh et al. (|2C)0C| ) 
Henry et al. (|2000^ 



Mayor & Q ueloz ( |1995| ) 
Vogt et al. (^000|) 



mined from Hipparcos data - ESA 1997 ), and within 20 pc 
with {B — V) < 0.9. This sample includes about 50 stars, 
most of which were observed. Some of the stars in this 
sample have known planetary systems, namely HD 1237, 
HD 13445, HD 17051, HD 22049, HD 28185 and HD 217107 
(HD 192263, would also have been included if we had not 
imposed any limit on {B — V)). All these objects were 
excluded from our final sample (43 objects) presented in 
Table |l} 

In what concerns the stars with planets, all known ob- 
jects having precise spectroscopic iron abundance determi- 
nations using a similar technique were used. This excludes 
G1876 (Delfosse et al |1998| ; Marcy et al. |1998D , HD 195019 
(Fischer et al. [1999 ) and the recently announced plan- 
ets ar ound HP 160691, H D 2744 2 and HD 179949 (Tinney 
et al. ^001| ; Butler et al. |200lD . Also, from the set of 9 
new planet host stars recently announced by the Geneva 



group (Mayor et al. 2001b) only two had available data 
(HD 28185 and HD 80606). The only exception was made 
respecting BD— 10 3166, which was include d in th e planet 
surveys for its high [Fe/H] (Gonzalez et al. 1999). All the 
other stars were part of search programmes for planets in 
volume-limited samples (Udry et al. 2000; Marcy et al. 



^oooD , and there is no particular reason to take out any of 
the targets (possible sources of bias are discussed below). 

2.2. Observations and data reduction 

For the "comparison sample" , spectra of S /N between 
150 and 350 were obtained using the CORALIE high- 
resolution spectrograph {R = AA/A ^ 50 000), at the 
Swiss 1.2-m Euler Swiss Telescope (La Silla, Chile) in 2000 
mid November. The spectra were reduced using IRAfQ 
tools, and the equivalent widths were measured by fitting 
a Gaussian function to each of the lines using the "k" key 
within splot. 

Some planet host stars were also observed with 
CORALIE (adding a few objects to the study presented 
in Paper I), some of them with no previous spectroscopic 
abundance determinations (see Table H) . Observations of 
some stars with planets were also complemented with 
spectra of S/N 300 taken with the FEROS spectrograph 
(i? - 48 000), at the ESO 1.52 m Telescope (La Silla), 



^ IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, U.S.A. 
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Fig. 1. Left: distribution of stars with planets (shaded histogram) compared with the distribution of field dwarfs 
presented in this paper (empty histogram). The dashed histogram represents the planet host sample if we consider 
only stars forming part of the CORALIE survey (see text for details). The vertical lines represent stars with brown 
dwarf candidate companions having minimum masses between 10 and 20Mjup. Right: The cumulative functions of 
both samples. A Kolmogorov-Smirnov test shows the probability of the stars' being part of the same sample is around 
10-7. 



during the nights of 2000 November 8 and 9. The UES 
spectrograph (i? ~ 55 000) at the 4 m WHT (La Palma, 
Canary Islands) was also used to obtain a spectrum for 
the planet host star HD 190228. 

2.3. Spectroscopic analysis 

The technique used has already been described extensively 
in Paper I. Abundances and atmospheric parameters were 
determined using a standard local thermodynamic equi- 
librium (LTE) analysis with a revised version of the line 
abundan ce co de MOOG (Sneden 1973 ), and a grid of 
Kurucz ( 1993 ) ATLAS9 atmospheres. As in Paper I, we 
used a set of iron lines taken from the list of Gonzalez 



and Fell lines are also tabulated. The masses were deter- 



& Laws (2000). Here, we added one more Fell line (the 
599lA presented in Gonzalez et al. 2001) and we excluded 
the previously used Fe ii line at 6432A fine, because it was 
giving systematically low values for most stars. 

The atmospheric parameters were obtained from the 
Fel and Fell lines by iterating until the correlation coef- 
ficients between loge(Fel) and xi, and between loge(Fei) 
and log {W\/X) were zero, and the mean abundance given 
by Fei and Fell lines were the same. This procedure gives 
very good results since the set of Fe i lines has a very wide 
range of excitation potentials. We used loge0(Fe) = 7.47. 

The complete results obtained for the comparison sam- 
ple are presented in Table |[ The errors in T^e, logp, 
and [Fe/H] for a typical measure are of the order of 50 K, 
0.15dex, O.lOdex, and O.OGdex respectively. The number 
of lines used for each star, and the dispersions for Fei 



mined from theoretical isochrones of Schaller et al. (1992), 
Schaerer et al. ( |1992| ) and Schaerer et al. (1993), using My 
computed from Hipparcos parallaxes and Tcs obtained 
from spectroscopy^. Only for HD 222237 we do not present 
a mass, given the high error in its determination. 

In Table ^ we present our results of the spectroscopic 
analysis for stars with low mass companions for which we 
have obtained spectra. Also included are the objects al- 
ready presented in Paper I and Santos et al. (2001), since 
the values were revised; the differences are always perfectly 
within the uncertainties (usually lower than 0.02 dex for 
[Fc/H]). Errors in the parameters were computed as de- 
scribed in Paper I. 

For a few stars in our sample we had both CORALIE 
and FEROS spectra. No important systematics are evi- 
dent in the results from the two spectrographs (due to 
errors in flatfields or background subtraction), and we sim- 
ply used a mean value of the results. 

We compared our [Fe/H] determinations for the eight 
stars presented in this paper in common with the stud- 
ies of Gonzalez et al . (|1999 ), Gonzalez & Laws ( ^000 ), 
and Gonzalez et al. (2001) to look for any possible sys- 
tematics. The stars and [Fc/H] obtained by these authors 
are HD1237 (0.16), HD 16141 (0.15), HD 17051 (0.19), 
HD 38529 (0.37), HD 52265 (0.27), HD 75289 (0.28), 
HD 210277 (0.24) and HD 217107 (0.36). The comparison 



^ The errors in mass represent simply formal errors, com- 
puted using the uncertainties in T^g and Mv. 
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Fig. 2. The shape of the distribution of metallicities for 
stars with planets increases with [Fe/H]. This represents 
also an increase with Z, the mass fraction of heavy ele- 
ments Z . 



shows that the formal mean difference in [Fe/H] obtained 
by these authors and from the present work is virtually 
0.00 dex; the rms of the differences is ^O.Oldex. The at- 
mospheric parameters for these stars are also usually the 
same within the errors. We thus conclude that we can in- 
clude their determinations for stars with planets in our 
analysis (as expected since our studies were similar in all 
aspects to theirs) without introducing any systematics, 
thus increasing the number of available planetary host 
candidates used in the current study. The list of stars with 
planets whose determinations were made by other authors 
is presented in Table 0. 



3. The metallicity distributions 

In Fig. 1^ we present the [Fe/H] distributions of both sam- 
ples described above. For the stars with planets we in- 
cluded both the objects presented in Tables || and |[ As 
is clear from the plot, stars with planets are significantly 
more metal rich than field stars without giant planet com- 
panions. While the stars-with-planets sample has a mean 
metallicity of -1-0.15 ± 0.23 the mean [Fe/H] for the field 
star sample is —0.10 ± 0.18 (here the errors represent 
the rms around the mean [Fe/H]). A Kolmogorov-Smirnov 
test (Fig. 1^, right) shows that the probability that both 
samples belong to the same population is about '^10~^. 

At this point it is important to discuss possible sources 
of bias. Is the star-with-planet sample completely un- 
biased? If we make the same plot using only the stars 



forming part of the CORALIE sample^ (from which the 
comparison stars were taken), the result is exactly the 
same (dotted histogram in Fig. |l]) : the same general shape 
and difference is found. In fact, this result cannot be re- 
lated to a selection bias since, as discussed above, the 
most important planet search programmes make use of 
volume-limited samples of stars. The only exception is 
BD-10 3166 (Butler et al. pOOOj), chosen for its high metal- 



licity. This star was not included in our analysis. It is 
worthy of note that the six planet-host stars that were in- 
cluded in our volume-limited sample (HD 1237, HD 13445, 
HD 17051, HD 22049, HD 28185 and HD 217107) have a 
mean [Fe/H] of -^0.12. 

No important systematics are expected concerning the 
magnitudes of the objects. On the one hand, for a given 
colour higher [Fe/H] stars are more luminous. Also, higher 
metallicity implies more and deeper lines, and thus a more 
precise determination of the velocity. But a star with more 
metals is also, for a given mass, cooler and fainter. For ex- 
ample, doubling the metallicity of the Sun (i.e. increasing 
[Fe/H] to ^ 0.30) would make its temperature decrease by 
more than 150 K, and its luminosity by a factor of ~1.2 



(Schaller et al. 1992; Schaerer et al. 1993). As we shall see 



below, for the mass intervals for which we have a good 
representation of both samples, stars with planets are al- 
ways significantly more metal-rich. Furthermore, at least 
in the CORALIE survey, exposure times are computed in 
order to have a photon-noise error at least as low as the 
instrumental errors. 

Given the uniformity of the study, we may conclude 
that the plot in Fig. |l| represents a proof that the stars now 
observed with giant planets are, on average, more metal- 
rich than field stars. For the record, it is also interesting to 
verify that the mean value of [Fc/H] obtained by Favata et 
al. (1997) for their volume-limited sample of G-dwarfs is 
—0.12, which means that there are almost no systematics 
between the method used by these authors and that used 
for the current study. 

Also remarkable in Fig. [| is the interesting position 
of the three low-mass (10Mjup<Msini<20Mjup) brown- 
dwarf candidates. Although it is too early to arrive at any 
firm conclusions, the position of one of them (HD202206), 
with [Fe/H] = -1-0.37 is strongly suggestive of a common 
origin with the lower-mass planets. On the opposite side 
of the distribution there is HD 114762 ([Fe/H] = -0.60), 
the most metal-poor object among all those studied in 
this paper. This "dispersion" might be interpreted as a 
sign that the frontier between brown dwarfs and massive 
giant planets is very tenuous (and probably overlaps) with 
regard to the mass limit; we are possibly looking at results 
of different formation processes (e.g. Boss 2000). 



These include HD 1237, HD6434, HD 13445 
HD 17051, HD 19994, HD 22049, HD 28185, 



HD 75289, 
HD 121504, 
HD 169830, 
HD 222582. 



HD 16141, 
HD 52265, 

HD 82943, HD 83443, HD 92788, HD 108147, 
HD 130322, HD 134987, HD 168443, HD 168746, 
HD 192263, HD 210277, HD 217107 and 
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Fig. 3. Panel (a): Metallicity distribution of stars with planets (dashed histogram) compared with the distribution 
of a volume limited sample of field dwarfs (empty histogram) - see text for more details; panel (b): correcting the 
distribution of stars with planets from the same distribution for stars in the volume-limited sample results in a 
even more steep rise of the planet host star distribution as a function of [Fe/H]. Note that the last bin in the planet 
distribution has no counterpart in the field star distribution; panel (c): same as (a) but for a field distribution computed 
using a calibration of the CORALIE cross-correlation dip to obtain the metallicity for ^1000 stars; panel (d): same 
as (b) when correcting from the field distribution presented in (c). 



4. How to explain it? 



ferences and to the implications this result might have for 
planetary system formation and evolution scenarios. 



The problem concerning the significance of the high- 
metallicity trend observed for stars with planets being 
"solved" , our enquiry is turned to the cause of these dif- 
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Fig. 4. Distribution for the field star sample (empty histogram) if we add 15 earth masses of iron (left) or a quantity 
of iron proportional to Z (right) - in the latter case, 0.003xZ. The distributions peak at about the same value as the 
distribution of stars with planets (dashed histogram), but have completely distinct forms and extents in [Fe/H]. See 
the text for more details. 



A.l. Comparing the distributions 

One very interesting feature of the distribution of stars 
with planets is well illustrated in Fig. |2[ The shape of 
this distribution is well described by a quadratic rising 
with [Fe/H] (up to [Fe/H] ^ +0.35 - dashed curve in the 
diagram). In other words, the distribution is rising with Z 
(the mass fraction of heavy elements) up to a value of 0.04, 
falling then abruptly. This sharp cut-off in the distribution 
of stars with planets is probably due to the fact that we 
may be looking at a limit on the metallicity of the stars 
in the solar neighbourhood. 

Clearly, the rise is expected if we consider that the 
probability of forming planets is proportional to the metal- 
licity content. In one sense, an increase in the number of 
dust particles (N) in the young disc will increase the prob- 
ability of a shock by iV^ . The rate of formation of planetes- 
imals, and thus the probability of creating a planet core 
with sufficient mass to accrete gas from the disc before it 
disappears, will increase quadratically with Z (considering 
that there is a linear relation between these two variables) . 

But this model is probably too simple. Other, more 
complicated mechanisms intervene, however, such as the 
fact that an increase in [Fe/H] will also probably increase 
the opacity of the disc, changing its physical conditions 
(e.g. Schmidt et al. 1997| ), such as the accretion rate or 
its vertical structure. The current result may be used in 
exactly this way in constraining these mechanisms. 

On the other hand, the distribution of [Fe/H] in a 
volume-limited sample of stars (with and without planet 
hosts) apparently decreases wi th in creasing [Fe/H] for 
[Fe/H] > 0.0 (e.g. Favata et al. |l997[). Taking this effect 



into account would result in an even steeper rise of the 
relative frequency of stars with planets with the metallic- 
ity. 

This is exactly what we can see from Fig. In pan- 
els (a) and (c) we present a comparison between the 
[Fe/H] distribution of stars with planets (shaded bars) and 
the same distribution for two volume-limited samples of 
dwarfs (empty bars). The volume- limited distribution in 
panel (a) was taken from this paper (objects presented 
in Table |l| plus the six stars with planets making part 
of our volume limited sample - see Sect. 2.1 for details). 
For panel (c), the distribution represents the metallicities 
of about 1000 stars in the CORALIE planet-search sam- 
ple as determined from a calibration of the CORALIE 
cross-correlation function surface (Santos et al., in prepa- 
ration)^. Since this two field star distributions represent 
the actual distribution of metallicities in the solar neigh- 
borhood (within statistical errors), we can use them to es- 
timate the "real" percentage of stars with planets for each 
metallicity bin that we would find if the metallicity distri- 
bution of field stars was "flat" - panels (b) and (d) . As we 



^ The cross-correlation dip ca n be seen as an av erage spectral 
line; as shown by Mayor (1980) and Pont (1997) its surface is 



well correlated with [Fe/H] and {B — V). A calibration of its 
surface with these two variables can in fact be used to obtain 
in a very simple way precise metallicity estimates, at least in a 
statistical sense. In this particular case, the calibration used to 
construct this "comparison" sainple uses as "calibrators" the 
stars presented in Tables |l| and y, and thus the results are in 
the same [Fe/H] scale as our spectroscopic values. The precision 
of the calibration is remarkable (rms~0.01 dex). The resulting 
distribution is quite symmetric, and peeks at ~— 0.05 dex. 
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can see from these plots, there is a steep rise of the percent- 
age of stars with planets with [Fe/H] (we will concentrate 
on panel (d) given that it presents better statistics than 
panel (b)). Actually, the plot shows that more than 75% of 
the stars with planets would have [Fe/H]>+0.2dex. But 
even more striking is the fact that ~85% of the surface of 
the distribution falls in the region of [Fe/H]>0.0. Although 
this is a preliminary result, the plots leave no doubts about 
the strong significance of the increase in the frequency of 
known planetary systems found with the metallicity con- 
tent of their host stars. 

We note that for the low metallicity tail of the dis- 
tribution, this result is based on relatively poor statistics. 
The small relative number of stars with planets found with 
low [Fe/H] values is, however, perfectely consistent with 
recent studies that failed to find any planetary host com- 
panions amid the stars in th e metal-poor globular cluster 
47Tuc (GiUiland et al. ^OOOl) . 



4.2. A simple "pollution" model 

We tried to determine whether it would be possible to ob- 
tain a metallicity distribution for the stars without planets 
similar to that of stars with extrasolar planets by simply 
adding material to their convective envelopes. To make 
this simple model, we computed the convective zone mass 
for each star using Eq. (5) and (6) of Murray et al. ( 2001 ) 
- metallicity-dependent relations for M < 1.2Mq; only 
stars in this mass range were used, since the comparison 
sample is not complete for higher masses. We are suppos- 
ing that the "pollution", if significant, has taken place af- 
ter the star reaches the ZAMS. As discussed in Paper I, if 
we consider the pollution to occur in a pre-main-sequence 
phase, higher quantities of material would be needed, since 
the convection envelope masses are bigger (D'Antona & 
Mazzitelli 1994). In fact, recent results show that gas-disc 
lifetimes can be higher than previously thought (up to 20 
Myr - Thi et al. pOOlj ), and thus we can probably expect 
significant pollution when the star is already on the main 
sequence. 

Two models were then constructed: the first was built 
on the supposition that the quantity of material falling 
into the star is the same for all the objects. A second, and 
possibly more realistic model, considered the pollution to 
be proportional to Z. The quantity of "added material" 
was chosen so that the distributions (stars with and with- 
out planets) peak at about the same value of [Fe/H]. 

The results are shown in Fig. ^ (open bars), compared 
with the results obtained for stars with planets within the 
same mass regime. As we can see, the shape of the planet 
sample, particularly the sharp cutoff at high metallicity, is 
not correctly obtained for any of the models. The result- 
ing samples are more symmetric, but most of all, in both 
models there are always a few points that end up with 
extremely high metallicities (> 0.8 dex). This is particu- 
larly unpleasant, since the inclusion of higher-mass stars 
(such as some of the observed exoplanet hosts) would even 




Fig. 5. A plot of [Fe/H] as a function of mass for planet 
host stars (filled circles) and "comparison" stars (open 
squares). The line represents the upper limit for the field 
star [Fe/H]. See text for more details. 



boost some bins to higher values, given that the mass of 
the convection zo ne is s upposed to be lower. As discussed 
by Murray et al. (2001), however, Li and Be observations 
suggest that the total "mixing" mass increases after 1.2 
Mq , eventually reducing this problem. Note also that the 
strong cut-off in [Fe/H] for the planet host star sample 
remains, even though this time we have not included stars 
with M > 1.2 Mq. This shape is not at all reproduced 
by our computations. 

We caution, however, that we are using very simple 
models to make this calculation. The inclusion of more 
parameters, such as a time- or disc-mass-dependent ac- 
cretion (which might be related to the stellar mass), or a 
ma ss-de pendent size of the central disc hole (e.g. Udry et 
al. |200l| ), might modify the results. In any case, we have 
shown that a simple "pollution" model cannot explain the 
observations. 



4.3. The mass dependence 



It has been suggested (Laughlin & Adams 1997) that if 
the pollution scenario is correct one would be able to see 
a trend of [Fe/H] with stellar mass, since the high-mass 
stars have tiny convective envelopes, and thus a higher 
probability of having a noticeable increase in [Fe/H] (here 
we have chosen stellar mass and not the mass of the con- 
vective shell, since if the pollution scenario is correct, the 
original mass for the convection zone for the stars with 
planets would not be correctly calculated). In Fig. ^ we 
present such a plot. 
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There are three features in the plot that deserve a 
comment. First, the region for M > 1.2 Mq has a very 
fow number of stars without planets. This is because 
in a volume- limited sample of stars, there are very few 
dwarfs in this mass regime (the IMF decreases with mass) . 
However, these stars are bright, and thus easier targets for 
planet searches, which explains the high number of plan- 
etary candidates in this region. On the other hand, for 
masses M < 0.8 Mq, we see the opposite effect: stars 
are fainter, and thus difficult targets for planet searches, 
but are more numerous in a volume-limited sample. There 
are probably, however, two more biases in the diagram. 
In fact, the comparison sample was cut in {B — V). For 
a given temperature, a highcr-metallicity star also has 
higher mass. Thus, stars in the upper part of the dia- 
gram (more metal-rich) are moved to the righthand side 
relative to stars in the lower part. For a fixed Tcs, 0.2 
dex in metallicity imply ~ 0.05 AfQ. On the other hand, 
if we take a given volume-limited sample of dwarfs, lower- 
mass dwarfs tend to be older than high-mass dwarfs; in 
proportion there are thus fewer metal-rich dwarfs in the 
"high" -mass region of the plot. These facts might explain 
the lack of stars with planetary companions presenting 
M > 1.2M0 and low [Fe/H]. 

In any case, the most striking thing in the diagram is 
the fact that all comparison stars have [Fc/H] < 0.17 dex, 
while about 55% of stars with planets have metallicity 
above that value. Furthermore, contrary to former studies 



Table 4. Planetary and orbital parameters used in Fig. |^ 
and 0. Stars with multiple companions are not included. 
The values were compiled from the literature. 



(Laughlin 2000) this diagram does not suggest that the 
metallicity of the upper envelope of stars with planets is 
increasing with stellar mass. This point renders unlikely 
the possibility that the high metallicity of these stars is 
due to pollution. As discussed above, however, for stellar 
masses higher than ~ 1.2 Mq, the total "mixing" mass 
may increase. This may explain to some extent the lower 
value of the upper envelope of [Fe/H] for masses above 
this limit, but it is not certain that it is enough to explain 
the difference. An extension of our comparison sample for 
stars with M > 1.2 Mq would be very useful. 

As mentioned, our results and interpretations are op- 
posed to those presented by Laughlin (200C), who used 
photometric indices to compute the metallicities for a sam- 
ple of stars with planets, and for a large volume-limited 
sample of dwarfs used as "comparison" . We believe that 
there are two important reasons for the difference. First, 
the trend he finds in [Fe/H] vs. stellar mass is mainly due 
to the few "low" -[Fe/H] planet host stars having M < 1 
Mq. On the other hand, in his plot there are a few very 
high metallicity "comparison" stars in the region of M < 
1.1 Mq. If we believe in the results presented in this pa- 
per, the probability that these stars harbour planets is 
very high. His sample might in fact be biased by the pres- 
ence of planet host stars in his "comparison" sample. 

If we concentrate on the mass interval for which we 
have both stars with planets and comparison stars (region 
with 1.2 Mq > M > 0.7 Mq), we obtain a mean [Fe/H] = 
-1-0.13 for the star with planet sample, and —0.08 for the 
comparison star sample. This makes a difference of about 
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4.23 


0.05 


0.00 


217107 


HD 217107 


1.34 


7.11 


0.071 


0.14 


222582 


HD 222582 


5.44 


575.9 


1.35 


0.71 



0.21 dex, for a mean stellar mass of 0.95 M0. To pollute 
the convective envelope of a 0.95 solar mass star contain- 
ing about 0.04 Mq in the convection zone (D'Antona & 
Mazzitelli 1994) in order to obtain the observed differ- 
ence, one would need about 10 earth masses of pure iron. 
Considering the composition of CI chondrites, this would 
mean ~ 5 times more in silicate material, a value that 
seems excessively high. The addition of giant planets, also 
very rich in gas, would need even more material. For exam- 
ple, the addition of Jupiter to the present-day Sun would 
only be able to increase its metallicity by ~ 0.05 dex, while 
adding two jupiters would increase 0.08 dex (e.g. Gonzalez 
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Fig. 6. Left: distribution of [Fc/H] for stars with planets orbiting with semi-major axes lower than 0.1 AU (shaded 
histogram), and with semi-major axes greater than 0.1 AU (open histogram). The vertical lines denote stars with more 
than one planetary-mass companion. Right: cumulative fraction of the [Fe/H] distribution for the long (thin line) and 
short (thick line) period systems. Althought some trend is suggested regarding a possible higher metallicity for the 
stars with short period planets, a Kolmogorov-Smirnov test gives a probability of ^^0.7, which is not significant. 



1998). We are supposing here that Jupiter has a rocky 
core, which is not necessarily true (Guillot |1999 ). 



Murray et al. (2001) have reported evidence that all 
stars in the solar neighbourhood have accreted about 
0.4 Earth masses of iron, suggesting that "terrestrial- 
type material is common around solar type stars." This 
seems to be particularly evident in the objects with M > 
1.3-1.4Mq, which is in accordance with expectations 
(Laughlin et al. 1997 ). If confirmed, and as discussed by 
these authors, this result further stresses the difficulty of 
obtaining the [Fe/H] differences we observe considering a 
pollution model: 0.4 Earth masses of iron are definitely 
not able to explain the observed differences between stars 
with planets and stars "without" planets. 



4.4. Using evolved stars as a test 

One other interesting note can be added by looking at the 
"star-with-planet" sample. If the pollution scenario were 
correct, we would expect evolved stars (sub-giants) hav- 
ing planetary companions to have lower metallicities than 
their dwarf counterparts, since the convection envelope in- 
creases in mass as the star evolves off the main sequence 



(Sackmann et al. 1993), diluting its metallicity content. 

There are eight stars in Tables || and || having spec- 
troscopic loggr values below ^ 4.10 dex (probably already 
slightly evolving away from the main sequence): HD 10697 
(3.96), HD 38529 (4.01), HD 117176 (3.90), HD 143761 
(4.10), HD 168443 (4.10), HD 169830 (4.04), HD 177830 



(3.32) and HD 190228 (4.02). If we compute the mean 
[Fe/H] we obtain a value of -1-0.08 ± 0.25 (with values 
going from —0.29 to +0.39), only slightly lower than the 
value of -1-0.15 ± 0.22 for the sample comprising only stars 
with logg > 4.10 (values going from —0.55 to -1-0.50). As 



mentioned also by Gonzalez et al. (2001), this is even more 
striking when we notice that the higher-metallicity star in 
this "low logg" sample (HD 177830) is that having the 
lower surface gravity. 

5. Metallicity and orbital parameters 

Given that we already have more than 40 extrasolar plan- 
ets with high-precision metallicity determinations, we can 
start to think about looking for possible trends in [Fe/H] 
with planetary mass, semi-major axis or period, and ec- 
centricity. 

Gonzalez (1998) and Queloz et al. ( ^000 ) have sug- 
gested that stars with short-period planets (i.e. small 
semi-major axes) may be particularly metal-rich, even 
amongst the planetary hosts. The number of planets that 
were known by that time was, however, not enough to 
arrive at a definitive conclusion. Here we revisit that dis- 
cussion with much better statistics. 

In Fig. ^ we present a plot of the distribution of planets 
for two orbital "radius" regimes: stars with planets orbit- 
ing with semi-major axes lower than 0.1 AU (shaded his- 
togram), and with semi-major axes greater than 0.1 A.U 
(open histogram). The values were taken from Table |[ 
Stars with multiple companions are denoted by the verti- 
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cal lines and were not included in the histograms. These 
include, among the stars for which we have [Fe/H] values 
available, the multiple systems around v And (Butler et al. 



199£), HP 83443 (Mayor et al. 2001a), HP 168443 (Udry 



et al. pOOlD , and HP 82943 (Mayor et al. |2001b|). 

Some trend can be seen in the plot (right panel) sug- 
gesting that stars with short period planets might be 
more metal-rich then their long-period counterparts. A 
Kolmogorov-Smirnov test, however, gives a probability 
that both distributions make up part of the same sample 
of '^75%. We have tried to make the plot using different 
limits of the semi-major axes; no further conclusions can 
be drawn. It is interesting to note that one of the longest- 
period planets discovered up to now (around the star 14 
Her) has a value of [Fe/H] = +0.50. 

Similar conclusions are drawn respecting the eccentric- 
ity or the companion mass (Fig. Q). The current results do 
not suggest any clear difference in [Fe/H] between high- 
and low-mass/eccentricity planets. If there is some trend, 
then more data (and probably even more precise abun- 
dances) are needed to permit conclusions. In the same 
way, no conclusions are possible regarding the position of 
the planetary systems in the plot (Fig. ||). 

6. Concluding remarks 

We have presented for the first time a uniform and un- 
biased comparison between the metallicity distributions 
of a sample of stars with planets and a sample of stars 
"without" planets. The main results can be summarized 
as follows: 

— The currently known stars with planets are substan- 
tially metal-rich when compared with non-planetary 
host dwarfs. The mean difference in [Fe/H] is '--^0.25 dex 
and is clearly significant. 

— The shape of the metallicity distribution of stars with 
planets has a very clear rise with [Fe/H], i.e. with Z, 
the mass fraction of heavy elements. This result, al- 
though still based upon "only" 44 objects, may help 
put constraints on the models and conditions leading 
to giant-planet formation. 

— The clear cut-off seen for the star-with-planet distribu- 
tion at [Fe/H] > -1-0.5 may be interpreted as a "limit" 
for the metallicity of solar-type dwarf stars in the so- 
lar neighbourhood; it is very difficult to explain it in 
terms of a convective envelope "pollution" scenario. 

— The metallicity of stars with planets cannot be ex- 
plained by a simple "pollution" model. The differ- 
ences in [Fe/H] between stars with planets and field 
stars does not seem to be explained by such mecha- 
nisms, but needs an original (cloud) metallicity "ex- 
cess" . Otherwise huge quantities of iron-rich material 
would have to be engulfed, which is not easy to explain 
in face of the current observations. Our results do not 
rule out, however, planet accretion as an ongoing pro- 
cess in the early stages of planetary formation (during 
which the convective envelopes are more massive). 



— An analysis of the planetary orbital parameters does 
not reveal any clear trends with [Fe/H]. This point 
might give further support to the "primordial" ori- 
gin of the high metallicity of stars with planets (e.g. 
Gonzalez |1998| ; Pel Popolo et al. |2001| ). On the other 
hand, we cannot exclude the possibility that some 
trend may appear in the future, e.g. when some real 
solar system analogues are found. 

The current result concerning the high metallicity of 
stars with planets may represent the simple fact that the 
higher the metallicity of the star, the higher will be the 
probability that planet formation will occur. 



Recently, Israelian et al. (2001) have discovered that 
the planet host star HP 82943 has ^Li in its atmosphere. 
This discovery is interpreted as evidence of the infall of 
a planet into the star, most probably by planet disrup- 
tion (Rasio & Ford |1996D . As discussed by the authors, 
however, the high [Fe/H] value of this star is not neces- 
sarily related to the infall of a planet, as this would not 
be able to enhance the [Fe/H] value by more than a few 
hundredths of a dex. This is exactly what w as fou nd for 
the pair 16Cyg A and B (Laws & Gonzalez (2001). 

In fact, the present results do not exclude the possibil- 
ity that pollution may play a role (eventually important 
in some cases, such as for the most massive dwarfs), but 
rather that it is not the key process leading to the observed 
high metallicity of the planet host stars. Recent work by 



Murray et al. (2001) supports the idea that pollution, on 
a small scale, is very common among the stars of the so- 
lar neighbourhood. If confirmed, however, this result does 
not change in any way the conclusions of this paper, since 
their models predict "pollution" of the order of 0.4 Earth 
masses of iron, only noticeable in stars with M >1.5Mq. 

One remarkable point that becames evident from the 
current work is the fact that our Sun occupies a "modest" 
position in the low [Fe/H] tail of the metallicity distribu- 
tion of stars with planets. A look at the orbital parame- 
ters in Table |^, however, tells us that, to date, no "real" 
Solar System analogs were found. This lead us to specu- 
late about possible different formation histories for these 
systems. Note, however, that we cannot draw any conclu- 
sions until other Solar System analogs are found. 

The future addition of more stars in both samples stud- 
ied here, and the detailed analysis of other chemical ele- 
ments (in particular the light elements Li, Be and B) are 
of particular interest for the future, both to consolidate 
current results, and to understand better the multitude of 
different planetary systems known today. The direct com- 
parison of other elemental abundances (e.g. a-elements, 
C, or O) between stars with planets and "single" field 
dwarfs will, on the other hand, not be an easy task given 
the apparent lack of stars with high [Fe/H] having no giant 
planet companions. Meanwhile, we expect that theoretical 
models of disc evolution and (giant) planetary formation 
will be constructed to explain in a detailed way the physics 
beyond the current results. 
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Fig. 7. Left: distribution of [Fc/H] for stars with planets having masses lower than 1 Afjup (shaded histogram), and 
greater than 1 Afjup (open histogram). Right: Similar to the left panel, but for eccentricity (limit at 0.1). In both cases, 
a Kolmogorov-Smirnov test results in a non-significative probability of being two different populations. 
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